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Forecast of prices of financial assets including gold is of considerable 
importance for planning the economy. For centuries, people have been 
holding gold for many important reasons such as smoothening inflation 
fluctuations, protection from an economic crisis, sound investment etc.. 
Forecasting of gold prices is therefore an ever important exercise undertaken 
both by individuals and groups. Various local, global, political, psychological 
and economic factors make such a forecast a complex problem. Data analysts 
have been increasingly applying Artificial Intelligence (AI) techniques to 
make such forecasts. In the present work an inter comparison of gold price 
forecasting in Indian market is first done by employing a few classical 
Artificial Neural Network (ANN) techniques, namely Gradient Descent 
Method (GDM), Resilient Backpropagation method (RP), Scaled Conjugate 
Gradient method (SCG), Levenberg-Marquardt method (LM), Bayesian 
Regularization method (BR), One Step Secant method (OSS) and Broyden- 
Fletcher-Goldfarb-Shanno Quasi Newton method (BFGS). Improvement in 
forecasting accuracy is achieved by proposing and developing a few 
modified GDM algorithms that incorporate different optimization functions 
by replacing the standard quadratic error function of classical GDM. Various 
optimization functions investigated in the present work are Mean median 
error function (MMD), Cauchy error function (CCY), Minkowski error 
function (MKW), Log cosh error function (LCH) and Negative logarithmic 
likelihood function (NLG). Modified algorithms incorporating these 
optimization functions are referred to here by GDM_MMD, GDM_CCY, 
GDM_KWK, GDM_LCH and GDM_NLG respectively. Gold price 
forecasting is then done by employing these algorithms and the results are 
analyzed. The results of our study suggest that the forecasting efficiency 
improves considerably on applying the modified methods proposed by us. 
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1. INTRODUCTION 


Gold price plays an important role in any economic system. Historically, gold has been considered a 
‘frontier-less currency’ that may be traded at any time and under virtually all circumstances. Gold has 
conventionally been the most effective way to collect cash during the world wide stock exchange crash in 
1987, and again in 1997 and 1998 during the Asian crises. The prices of gold and other financial assets like 
equities and other commodities are often closely correlated. For instance, the correlation between gold and 


Journal homepage: http://ijai.iaescore.com 


Int J Artif Intell ISSN: 2252-8938 0 47 


equities may generally be negative since during the stock market booms funds are frequently transferred from 
gold to stocks and vice versa at times of crises. The correlation between gold and oil is generally positive and 
international disputes may increase the prices of oil and gold simultaneously. Accurate gold price forecasting 
is therefore very important to have a prudent investment policy which also helps mitigating financial risks. 
However, there is no mathematical formulation or model available to describe movement of gold price. In 
such a scenario, methods based on application of AI are inherently capable of making suitable forecasts. The 
advantages of AI methods lie in their ability to model complex and nonlinear relationships without any prior 
knowledge of the processes generating them [1]. AI methods are traditionally classified in three categories: 
Artificial Neural Network (ANN), Fuzzy logic (FL), Genetic algorithm (GA) and their hybrids which are 
combinations of these methods. ANNs have become preferred tools for forecasting of the financial markets in 
general. This is because of their ability to extract trends and patterns in a given data set. Trained ANNs can 
be considered as expert systems in their domain [2]. However, within various ANNs, choosing the best one 
for a given problem is a very important task [3]. The best choice is naturally the one giving most efficient 
forecast for the problem at hand. 

A large amount of work has been done by researchers demonstrating the forecasting power of 
AI/ANN based methods in various disciplines. A brief review of it is given here. A crisp review of the 
efficiency of AI techniques for diagnostic systems [4]. Various ANN methods like GDM, Conjugate 
Gradient, Quasi Newton and LM method to predict and classify the patients with heart disease [5]. Employed 
a backpropagation method with modified learning rate and momentum factors and applied the method to 
breast cancer and iris data sets [6]. Their study shows an impressive decrease in the training times of 
backpropagation which is otherwise inherently plagued by slow learning and local minima problems. A feed 
forward ANN model with the application of LM learning algorithm for short-term load forecasting of daily 
peak load and demonstrated that the proposed ANN model gives more accurate predictions with optimal 
number of neurons in the hidden layer [7]. ANN and Random Forest (RF) approaches for short-term 
photovoltaic (PV) output current forecasting (STPCF) for the next 24-hours. The results have shown that 
both the proposed techniques are able to perform forecasting of future hourly PV output current efficiently 
[8]. A merged Long Short-term Memory (LSTM) for forecasting ground visibility at the airport by 
combining time series of predictor variable with another moderating variable. Forecasting accuracy thus 
improved over the traditional LSTM method [9]. ANNs to diagnose melanoma skin cancer at an early stage 
with a high degree of accuracy [10]. Applied various neurocomputing methods for rainfall forecasting giving 
highly accurate results [11]. Proposed and tested a new hybrid model, i.e. Generalized Space-Time 
Autoregressive with Exogenous Variable and a Neural Network (GSTARX-NN) model for forecasting space- 
time data with calendar variation effect. They found that the hybrid GSTARX-NN model gave more accurate 
forecasts than the traditional GSTARX models [12]. Compared a few variants of exponential and 
backpropagation ANN models to forecast rice production. The results showed that neural network method is 
preferable to the statistical method since it results in lower Mean Square Error (MSE) and Mean Absolute 
Percentage Error (MAPE) [13]. AI has also been used in the financial sector for various kinds of forecasts, 
analyses and decision making. Used of machine learning algorithms is explored to analyse effect of financial 
news on stock market prices. Support Vector Machine (SVM) and RF algorithms were used and it is 
concluded that the RF algorithm gives better accuracy in comparison with SVM algorithm [14]. Proposed a 
procedure to construct Triangular Fuzzy Number from single point data and then used an Autoregressive 
model to forecast currency exchange rates of Association of South East Asian Nation (ASEAN) countries 
[15]. A Fuzzy Neural System (FNS) to forecast the inflation rate obtaining better results in terms of RMSE 
than traditional methods [16]. A fuzzy sets method predicting the Russia Trading System (RTS) index with a 
confidence level of about 90% though having only an incomplete data set at hand [17]. A hybrid model, 
wavelet radial bases function neural networks (WRBFNN), for forecasting of non-stationary time series and 
found that, in terms of MAPE and MSE, the proposed model is superior to the traditional wavelet feed 
forward neural networks (WFEN) model [18]. 

In the present work forecast of gold prices in the Indian Multi Commodity Exchange market, MCX, 
has been taken up. A comparative analysis of gold with other investment options helps one to prioritize 
capital investments [19]. Traditional mathematical models such as Moving averages (MA), Auto-Regressive 
Moving Averages (ARMA), Auto-Regressive Integrated Moving Averages (ARIMA) etc. as well as methods 
based on AI such as the ANN, Adaptive Neuro Fuzzy Inference Systems (ANFIS) etc. have been deployed 
by many researchers, the world around, for gold price forecasting achieving varying degree of success. After 
making a brief review of these results we present our work on gold price forecasting. 

In our work we have proposed a few novel methods based on the classical GDM approach belonging 
to the ANN class and developed algorithms based on them on MATLAB platform. For this purpose the 
classical optimization function, the quadratic error function, of the classical GDM method is replaced in turn 
by several other optimization functions. These newly developed algorithms, along with many other popular 
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ANN algorithms, are then deployed to forecast gold prices in India. This is done by systematically deploying 
all the above mentioned methods on the same data set of gold prices and computing various efficacy 
parameters of forecast including the most significant one the mean percent accuracy (MPA). An inter- 
comparison of the results obtained with various methods shows that the hereby proposed methods result in 
significant improvement in price forecasts. The most efficient optimization function is thus identified for the 
given data set and the ANN method employed. The manuscript is organized as follows: Section 2 presents a 
brief description of the research methodology adopted. Section 3 discusses the results obtained and finally 
Section 4 concludes the findings of the present work. 


2. RESEARCH METHODOLOGY 

A few ANN methods, classical and proposed ones, have been employed here in the present work. A 
brief description of these methods follows. The backpropagation (BP) algorithm that was introduced by 
Rumelhart is a well-known method for training multilayer feed-forward artificial neural networks 
[20]. It makes use of the classical GDM technique. The Multi-Layer Perceptron (MLP) networks 
trained by BP algorithm have been a popular choice in financial applications [21]. The BP 
algorithm suffers from slow convergence and traps in local minima of error function. Several 
modifications have been introduced by researchers for improving the local minima problem. These 
are: gradient descent with momentum, gradient descent with adaptive learning rate, gradient 
descent with both momentum and adaptive learning rate and resilient backpropagation [22]. Some 
other methods are based on variations of numerical optimization technique itself resulting in 
advanced ANN methods like the LM, Quasi Newton and Conjugate Gradient methods. Conjugate 
Gradient algorithms perform searches along conjugate directions requiring calculation of second 
derivatives. Nevertheless, its convergence is generally faster than that of the steepest descent 
direction method [23]. Quasi Newton method may converge still faster than conjugate gradient 
methods since it avoids calculation of second derivatives. It approximates the Hessian matrix and 
updates it at each iteration. The Broyden-Fletcher-Goldfarb-Shanno (BFGS) method has been one of the 
most successful methods in the class of Quasi Newton methods [24-25]. The one step secant (OSS) is 
formulated in such a way that it attempts to overcome the limitations of the techniques like Quasi Newton 
and conjugate gradient methods. The method has an inherent computational advantage since Hessian matrix 
storage at every iteration is not required [26]. The LM method interpolates between the Gauss-Newton 
and the gradient descent methods. It also avoids calculation of the Hessian Matrix by replacing it 
by the Jacobian needing lesser computation [27]. Some details of these methods can be found in the 
work done by [28]. Bayesian regularization minimizes a linear combination of squared errors and 
weights and it exhibits good generalization qualities [29]. 

In all the methods mentioned above the optimization function to be minimized, E, is the 
quadratic sum of the errors, the difference between the actual (ya) and forecasted (y») values of the 
variable of interest and is given by (1). 


2 
E=05 Yhilva-yp), = 0-5 0% er? (1) 


It is envisaged in the present work that variants of the objective function E may improve 
the forecasting efficiency. A total of seven classical methods mentioned earlier and five modified 
GDM algorithms proposed and developed by us, using different optimization functions, were 
chosen for the present study. Various optimization functions are discussed in next section. 


2.1. Modified GDM methods 
In the present work, five modified GDM methods/algorithms have been proposed and developed by 
us which incorporate different optimization functions instead of the standard quadratic error function in the 
classical GDM to forecast the gold price [30-31]. These five variants of the optimization function are 
discussed. 
1. Mean median (MMD) error function: Mean Median error function is given in (2). This function has the 
advantage of both the Mean error function and Median error function. It reduces the influence of large 
errors. 


E= Y,(2*(J1+ (e2/2)) - 1) (2) 


2. Minkowski (MKW) error function: Minkowski error function is given in (3). MSE error can be 
visualized as a special case of it with r =2. Variations of r around 2 have been observed to have 
significant influence on accuracy of forecast. Here r has been chosen to be 0.4. 
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E=Ynleil” (3) 


3. | Logcosh (LCH) error function: 
Log cosh error function is given in (4). It approximates x?/2 and |x| for small and large values of x 
respectively. 


E = Y,, In(cosh( e;”)) (4) 


4. Cauchy (CCY) error function: 
Cauchy error function is given in (5). Cauchy error function is known to show robustness against outliers. 


E=Y, c?/2In(4 + (2)") where c = 2.38 (5) 


5. | Negative Logarithmic Likelihood (NLG) error function: 
Negative Logarithmic Likelihood error function is given in (6) Negative logarithmic likelihood error function 
is a measure of accuracy of a classifier and is defined as follows. 


E = —=Y,,In(¢;) (6) 


The modified algorithms incorporating the above functions are here referred to as GDM_MMD, 
GDM_MKW, GDM_LCH, GDM_CCY and GDM_NLG respectively. Forecasting efficiency is judged here 
by examining various efficacy parameters described briefly next. 


2.2. Forecasting efficacy parameters 

Many efficacy parameters can be defined to measure degree of success in forecasting. Let ya denote 
the actual stock price and y; the forecasted one. The error in forecast is then given by e= (ya — yr). Let ‘n’ 
denote the total number of input-output sets, constructed from the time series of stock prices. Let the mean of 
all the actual stock prices being forecasted (the desired outputs) be denoted by Ymean. The set of efficacy 
parameters of forecasting considered here to investigate efficiency of AI methods are discussed [32]. 
1. Mean Error (ME) 


1 n 
ME=!¥",¢, (7) 


Mean error (ME) is a measure of the average difference between the actual (desired) prices and the 
forecasted ones and is given in (7). As it shows the direction of error it may also be called forecast bias. ‘ME’ 
may turn out to be deceptively small as positive and negative values of individual e, may cancel each other. 


2. Mean Absolute Deviation (MAD) 
1 n 
MAD = ~Y el (8) 


This parameter is sum of absolute values of errors and is given in (8). Positive and negative errors no more 
cancel each other. It does not however indicate the direction of errors. For a good forecast ‘MAD’ should be 
made as small as possible. 

3. Mean Square Error (MSE) 


MSE="52 (e) (9) 


It is clear from the above expression that ‘MSE’ is the average of squares of forecast errors and is given in 
(9) Large individual errors affect it significantly. This parameter does not give any idea of the direction of the 
overall error. 

4. Root Mean Square Error (RMSE) 


RMSE =/MSE (10) 


RMSE contains all the properties of MSE and is given in (10). It is widely computed and reported in 
literature. 
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5. Mean Percent Forecast Error (MPE) 


MPE=~Y", (=) *100 (11) 

n dy 
MPE is independent of the magnitude of the prices as it is the average of percentage of forecasting errors and 
is defined in (11). It indicates the direction of error. Individual errors of opposite signs nullify each other 
making it deceptively small. Like ‘ME’ its value should also be kept as small as possible for a good forecast. 
6. Mean Absolute Percent Forecast Error (MAPE) 


MAPE =-yn, | (=) *100| (12) 
t 


MAPE gives the average of percentage of absolute values of individual errors in forecast and is given in (12). 
It is independent of the scale of the data. Individual errors of opposite signs do not cancel each other. 
Sufficiently small value of ‘MAPE’ may ensure a good forecast. 

7. Mean Percent Accuracy (MPA) 


MPA = 100 — MAPE (13) 


‘MPA’ depends on ‘MAPE’ directly. It is defined in (13). As it gives mean percent accuracy and not the 
error, it may sound more pleasing to some users. 
8. Coefficient of Determination (R?) 


ERi(va-va). 
2 


TE (va,- Yeisen) 


R15 (14) 


R’ is a measure that represents the goodness of fit of a model and is given in (14). Large values of R? may 
suggest a good fit to historical data. It may however not be so while forecasting out of sample data [32]. 


2.3. Data collection and preprocessing 

Weekly time horizon has been considered in the present work. The gold price data is downloaded 
from a reliable site on investing [33]. A data window of 4 years of gold prices, from Jan 2015 to Dec 2018, 
has been considered here. The data consists of the prices of gold for the first trading day of every week in the 
entire time period. Prices are normalised between 0.1 and 0.9. Time series formation is done by making sets 
of weekly prices of 5 weeks as input and the price of the 6" week as the corresponding output. Time series so 
formed are divided in three categories: training, testing and validation. Training and testing series constitute 
about initial 80% of the total. About 70% of these series are randomly selected for training the network and 
the remaining 30% are used for testing. Validation series correspond to chronologically the latest data and 
constitute about 20% of the total.The validation data is thus not contained in the training and testing data and 
is referred to, here, as out of sample data. Forecasting accuracy will be examined for varying validation 
periods. 


3. RESULTS AND DISCUSSION 
In the present study the one week ahead price of gold is forecasted from the previous 5 weeks’ 

prices using the classical ANN methods and the modified GDM methods for various ANN architectures. The 

work done is categorized as below and discussed in detail in sections 3.1 to 3.3. 

1. Forecast of weekly prices of gold using classical ANN algorithms namely GDM, RP, SCG, LM, BR, 
BFGS, and OSS. 

2. Implementation of the proposed modified GDM methods to increase forecasting efficiency and inter- 
comparison of the results obtained with those obtained by classical GDM. 

3. Comparison with the work done by other researchers for gold forecasting using traditional ANN 
methods. 


3.1. Forecasting using classical ANN algorithms 


Forecasting of weekly gold prices has been done using the classical ANN algorithms namely GDM, 
RP, SCG, LM, BR, BFGS, OSS with the architecture 05-03-01. Efficacy parameters of forecasting obtained 
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for testing and validation phases are summarized in Tables 1 and 2. Important observations are summarized 
below each table. 

Figure 1 plots MPA, MSE and R? obtained in the testing phase for all the algorithms. It can be 
observed from Table | and Figure | that, in the testing phase, the maximum MPA (90.35%) is obtained by 
deploying RP algorithm and the least MSE (2.68E-03) and the maximum R? (0.9318) are obtained 
by BR algorithm. 


Table 1. Efficacy parameters of forecasting of gold prices using ANN (05-03-01) in testing phase 


ANN Method ME MAD MSE RMSE MPE% | MAPE% MPA% R? 
GDM 0.02150 0.0664 0.00603 0.7770 1.65 16.51 83.49 0.8261 
RP 0.00319 0.0409 0.00279 0.0528 -0.04 9.65 90.35 0.9289 
SCG -0.00324 0.0467 0.00470 0.0686 -0.64 10.91 89.09 0.8645 
LM -0.00213 0.0521 0.00393 0.0627 -1.35 12.08 87.92 0.8997 
BFGS 0.00484 0.0427 0.00295 0.0430 0.13 10.06 89.94 0.9248 
BR 0.00336 0.0406 0.00268 0.0517 -0.27 9.81 90.19 0.9318 
OSS 0.00231 0.0428 0.00289 0.0538 -0.05 9.96 90.04 0.9263 


MPA% in Testing Phase R? in Testing Phase 


MSE in Testing Phase 


= 


a 
=| a 


8 
[=2} 
ANN Method 


Figure 1. Comparison of MPA, MSE and Rin the testing phase of gold price forecast for classical methods 


Further, the results were validated on an out of sample data (validation phase) also for the 
chronologically latest 46 weeks. All the efficacy parameters except R? are shown in Table 2. R? is not a 
reliable performance indicator for the out of sample data. 


Table 2. Efficacy parameters of forecasting of gold prices using ANN 05-03-01 in validation phase 


ANN Method ME MAD MSE RMSE MPE% MAPE% MPA% 
GDM 0.0376 0.0450 0.00320 0.0566 4.61 5.72 94.28 
RP 0.0156 0.0372 0.00233 0.0483 1.92 4.86 95.14 
SCG 0.0148 0.0485 0.00375 0.0613 1.96 6.29 93.71 
LM 0.0179 0.0426 0.00318 0.0564 2.30 5.74 94.26 
BFGS 0.0096 0.0400 0.00274 0.0524 1.09 5.24 94.76 
BR 0.0145 0.0346 0.00190 0.0436 1.71 4.50 95.50 
OSS 0.0185 0.0400 0.00276 0.0526 2.19 5.20 94.80 


Figure 2 plots MPA and MSE obtained in the validation phase for all the algorithms. It can be seen 
from Table 2 above and Figure 2 that in the validation phase, the BR algorithm gives the best MPA (95.5%) 
and best MSE (1.90E-03). 
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MPA% in Validation phase MSE in Validation phase 


LM BFGS BR OSS RP SCG LM BFGS- BR 
ANN Method ANN Method 


Figure 2. Plots of MPA and MSE in the validation phase of gold price forecast for classical methods 


As the validation results are more significant than the test results for forecasting applications, BR 
can be considered as the best model for the data in consideration. Graphs in Figure 3, plot desired versus 
forecasted normalized prices of the testing phase and validation phase for BR algorithm with the ANN 
configuration 05-03-01. 


Testing Phase Validation Phase 


Event No. Event No. 
=== Desired 0/p (y) ==m@== Predicted o/p =="== Desired 0/p(y) | “====Predicted o/p 


Figure 3. Desired v/s forecasted output for gold price test data and validation data using BR with 
ANN 05-03-01 


3.2. Forecasting using the proposed GDM algorithms and inter-comparison with classical GDM 

Forecasting has also been done using the 5 novel modified GDM algorithms, GDM_MMD, 
GDM_MKW, GDM_LCH, GDM_CCY and GDM_NLG, discussed in section 2.1 with various architectures. 
These results are compared with those of the classical GDM results given in Section 3.1. Efficacy parameters 
obtained for testing and validation phases are summarized in Tables 3 and 4. 


Table 3. Forecasting performance during testing phase of classical GDM and modified methods of GDM 
with various configurations of ANN for gold 


ANN Method ANN Architecture MSE (Testing Phase) MPA% (Testing Phase) R? (Testing Phase) 
Classical 05-03-01 0.00600 83.49 0.8261 
GDM 05-05-01 0.00930 80.86 0.7314 
05-07-01 0.00610 86.58 0.8238 
GDM_MMD 05-03-01 0.00290 90.18 0.9162 
05-05-01 0.00330 89.80 0.9053 
05-07-01 0.00280 91.23 0.9192 
GDM_CCY 05-03-01 0.00398 87.14 0.8852 
05-05-01 0.00460 86.42 0.8668 
05-07-01 0.00440 89.14 0.8730 
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Table 3. Forecasting performance during testing phase of classical GDM and modified methods of 
GDM with various configurations of ANN for gold (Continue) 


ANN Method ANN Architecture MSE (Testing Phase) MPA% (Testing Phase) R’ (Testing Phase) 
GDM_KWK 05-03-01 0.00350 89.01 0.8986 
05-05-01 0.00300 80.42 0.9141 
05-07-01 0.00310 90.28 0.9092 
GDM_LCH 05-03-01 0.00400 88.41 0.8845 
05-05-01 0.00320 92.40 0.9078 
05-07-01 0.00280 90.76 0.9198 
GDM_NLG 05-03-01 0.00318 89.90 0.9189 
05-05-01 0.00572 88.01 0.8543 
05-07-01 0.00303 89.31 0.9229 


It can be observed from the Table 3 that in the testing phase the modified GDM methods give, in 
general, better accuracy than the classical GDM method for all configurations of the network. A comparison 
of all the MPAs obtained for various ANN architectures by using different algorithms, is shown graphically 
in Figure 4. 


MPA% in Testing Phase MSE in Testing Phase R? in Testing Phase 
0,010 


0,008 
0,006 
0,004 
0,002 
0,000 


05/03/01 05/05/01 05/07/01 05/03/01 05/05/01 05/07/01 05/03/01 05/05/01 05/07/01 
ANN Architecture ANN Architecture ANN Architecture 


=GDM =GDM_MMD mGDM_CCY =GDM mGDM_MMD BGDM_CCY | aco =GDM_MMD =mGDM_cCY 
=mGDM_KWK =GDM_LCH m=®GDM_NLG mGDM_KWK =GDM_LCH ==GDM_NLG | aco KWK mGDM_LCH mGDM_NLG 


Figure 4. Comparison of MPA, MSE and R? in the testing phase of gold price forecast for GDM and 
modified methods 


It can be concluded from Figure 4, that in the testing phase, the best MPA obtained is given by the 
GDM_LCH method with the architecture 05-05-01 and it equals 92.4%. The least MSE obtained is 2.80E-03 
and is given by GDM_MMD and GDM_LCH with the architecture 05-07-01. The maximum R? obtained is 
given by GDM_NLG with the architecture 05-07-01 and it equals 0.9229. Once the network is trained and 
tested it is expected to make good forecasts during the validation period. An important question can be raised 
here. How does the forecasting accuracy compare for various GDM methods, the classical and the modified 
ones and also how does it change as the validation time period is changed? In order to answer this question 
forecasts are made with all the GDM methods for varying lengths of validation period, namely, 4 weeks, 12 
weeks, 24 weeks and 46 weeks. The results are tabulated in Table 4. 


Table 4. Forecasting Performance during Validation Phase of classical GDM and modified methods of GDM 
with various configurations of ANN for gold 


Method ANN Forecasting Accuracy (Validation phase) 
Architecture 46 weeks 24 weeks 12 weeks 4 weeks 
MSE MPA% MSE MPA% MSE MPA% MSE MPA% 
Classical GDM 05-03-01 0.00320 94.28 0.00240 94.62 0.00170 95.39 0.00180 94.55 
05-05-01 0.01240 88.20 0.01100 87.89 0.00520 91.21 0.00300 93.85 
05-07-01 0.00940 89.78 0.00830 90.49 0.00890 88.84 0.00750 89.16 
GDM_MMD 05-03-01 0.00220 95.07 0.00150 96.08 0.00180 95.78 0.00360 92.97 
05-05-01 0.00270 94.58 0.00160 95.86 0.00150 95.83 0.00250 94.61 
05-07-01 0.00280 94.37 0.00230 94.68 0.00260 94.47 0.00450 92.23 
GDM_CCY 05-03-01 0.00156 95.94 0.00111 96.53 0.00102 96.36 0.00100 96.05 
05-05-01 0.00200 95.36 0.00160 95.79 0.00170 95.35 0.00200 94.01 
05-07-01 0.00380 93.26 0.00270 94.33 0.00410 92.53 0.00710 88.72 
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Table 4. Forecasting Performance during Validation Phase of classical GDM and modified methods of GDM 
with various configurations of ANN for gold (Continue) 


Method ANN Forecasting Accuracy (Validation phase) 
Architecture 46 weeks 24 weeks 12 weeks 4 weeks 

MSE MPA% MSE MPA% MSE MPA% MSE MPA% 

GDM_KWK 05-03-01 0.00250 94.84 0.00240 94.69 0.00180 95.49 0.00200 95.22 
05-05-01 0.00330 93.77 0.00300 93.98 0.00340 93.20 0.00350 92.24 

05-07-01 0.00260 94.65 0.00250 94.60 0.00300 94.16 0.00420 90.93 

GDM_LCH 05-03-01 0.00230 94.86 0.00160 95.96 0.00220 94.74 0.00430 91.94 
05-05-01 0.00270 94.58 0.00163 95.86 0.00150 95.83 0.00250 94.61 

05-07-01 0.00310 94.14 0.00240 94.72 0.00150 95.43 0.00210 94.43 

GDM_NLG 05-03-01 0.00230 94.96 0.00170 95.92 0.00150 95.96 0.00300 93.96 
05-05-01 0.01220 88.27 0.01110 87.99 0.00530 91.20 0.00340 92.69 

05-07-01 0.00223 95.23 0.00200 95.40 0.00180 95.92 0.00240 94.03 


It can be observed from the above table that all the modified GDM methods give better MPAs than 
the classical GDM method for all configurations of the network for the 46 weeks and 24 weeks validation 
periods. For 12 and 4 weeks validation periods also the modified methods give, in general, better MPA than 
the classical GDM. GDM_CCY outperforms all the other modified GDM algorithms as well as the classical 
GDM in terms of forecasting accuracy. It gives around 96% MPA across all validation time periods. For a 
better visualisation of our results a consolidated comparison of all the MPAs and MSEs obtained for various 
ANN architectures by using different algorithms for validation period of entire 46 weeks, is shown 
graphically in Figure 5. 
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Figure 5. Comparison of MPA and MSE in the validation phase of gold price forecast for classical GDM and 
modified methods 


It can be observed from the above figure that for this validation period, the GDM_CCY algorithm 
with the architecture 05-03-01, outperforms all other modified and classical GDM algorithms in terms of 
forecasting MPA and MSE. Graphs in Figure 6 plot the desired versus forecasted normalized prices in the 
testing phase and validation phase respectively for GDM and GDM_CCY algorithms with the ANN 
configuration 05-03-01. Table 5 shows all the forecasting efficacy parameters of these two methods for the 
testing phase. 


Table 5. Efficacy parameters of gold price forecasting using the classical GDM and GDM_CCY 
in testing phase 


Method ANN Architecture ME MAD MSE RMSE ~MPE% —_MAPE% —MPA% R? 
Classical GDM 05-03-01 0.0215 0.0664 0.00603 0.0777 1.65 16.51 83.49 0.8261 
GDM_CCY 05-03-01 0.0060 _0.0518 0.00398 _0.0631 -0.76 12.86 87.14 0.8852 
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Figure 6. Desired v/s Forecasted output for gold price test data and validation data using Classical GDM and 
GDM_CCY 


It can be also be observed that deployment of the modified algorithm, GDM_CCY, improves the 
MPA by ~ 4% over the classical GDM in the testing phase. Other parameters are also significantly better. 
Table 6 shows all the forecasting efficacy parameters (except R*) of the two methods for the validation 
period of 46 weeks. R? is not a reliable performance indicator for the out of sample data. 


Table 6. Efficacy Parameters of gold price forecasting using the classical GDM and GDM_CCY in 
Validation Period of 46 weeks 


Method ANN Architecture ME MAD MSE RMSE MPE% MAPE% MPA% 
Classical GDM 05-03-01 0.0376 0.0450 0.00320 0.0566 4.61 5.72 94.28 
GDM_CCY 05-03-01 0.0121 0.0312 0.00156 0.0394 1.40 4.06 95.94 


It can be seen from the above table that deployment of the modified algorithm, GDM_CCY, 
improves the MPA by about 1.6% over the classical GDM in the validation phase. Other parameters are also 
significantly better in the modified method. It can also be observed from Table 2 and Table 6 that in the 
validation phase, the modified algorithm GDM_CCY outperforms even the best performing classical 
algorithm i.e. the BR algorithm. 


3.3. Comparison with work undertaken by other researchers 

Various researchers have attempted to forecast gold price in the past usually employing the 
conventional ANN and other methods. Used recursive and rolling neural network models to forecast one- 
step-ahead sign variations in gold prices getting an average accuracy of 60.68% [34]. Forecasted gold prices 
with ANN and ARIMA and concluded that, based on three statistical criteria namely coefficient of 
determination (R’), mean absolute error (MAE), and root mean square error (RMSE), ANN outperformed 
ARIMA significantly [35]. Used ARIMA to forecast gold prices and obtained MAPE ranging between 3% 
and 5% [36]. Compared the gold price forecasting performance of an algorithm BAT-NN with that of 
ARIMA, ANN, ANFIS, MLP, RBF and GRNN. BAT-ANN produced smallest RMSE while ARIMA 
(1, 1, 3) gave the smallest value of MAPE 3.135 [37]. Observed that the PCA-GA-BP model gave better 
results than GA-BP and BP models, average relative errors of gold price prediction being 1.637%, 3.124% 
and 5.018% [38]. Considered monthly gold prices and compared forecasting performance of a classical ANN 
method with a hybrid model of ANN and GA. The hybrid model produced a lower MAPE (about 6.3) than 
that of the ANN model (about 7.6) [39]. 

In our study the algorithm GDM_CCY forecasts the weekly gold price with MPE about 96% 
(MAPE 4%) which compares well with the accuracies obtained by others summarized above. One should 
however note that a direct comparison cannot be done as the data set of the prices is not the same in all the 
above studies. The moot point here is that the modifications suggested by us result in significant 
improvements in forecasting efficiency. This is apparent when various efficacy parameters of forecast are 
compared with and without the proposed modifications. 
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4. CONCLUSION 

In the present work we proposed and developed a few novel ANN algorithms to enhance forecasting 
accuracy of gold prices in a popular Indian commodity market. For this purpose we proposed to employ 
several alternate optimization functions instead of the prevalent quadratic error function based on the mean 
square error. Weekly gold price forecasts were made using various classical ANN algorithms and the hereby 
developed modified ANN algorithms. Several efficacy parameters of forecast so obtained have thus been 
inter compared. It was thus demonstrated that significant improvement in several efficacy parameters 
including the most important one, the mean percent accuracy (MPA), can in general be achieved when the 
proposed alternate optimization functions were used. The inter comparison shows that the Cauchy error 
function based algorithm produced the best results among all the classical as well as other proposed 
algorithms for the data under consideration. It is therefore, hereby, suggested that the methods proposed by 
us like GDM_CCY and others should also be deployed before taking a final call on the price forecasts. 
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